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We present a theory of single-magnet flowmeter for liquid metals and compare it with experi- 
mental results. The flowmeter consists of a freely rotating permanent magnet, which is magnetized 
perpendicularly to the axle it is mounted on. When such a magnet is placed close to a tube carrying 
liquid metal flow, it rotates so that the driving torque due to the eddy currents induced by the flow is 
balanced by the braking torque induced by the rotation itself. The equilibrium rotation rate, which 
varies directly with the flow velocity and inversely with the distance between the magnet and the 
layer, is affected neither by the electrical conductivity of the metal nor by the magnet strength. We 
obtain simple analytical solutions for the force and torque on slowly moving and rotating magnets 
due to eddy currents in a layer of infinite horizontal extent. The predicted equilibrium rotation rates 
qualitatively agree with the magnet rotation rate measured on a liquid sodium flow in stainless steel 
duct. 

I. INTRODUCTION 

Flow rate measurements of liquid metals are required in various technological processes ranging from the cooling 
of nuclear reactors to the dosing and casting of molten metals. 1 Electromagnetic flowmeters are essential in the 
diagnostics and automatic control of such processes. A variety of electromagnetic flowmeters have been developed 
starting from the late 1940s and described by Shercliff. 2 The standard approach is to determine the flow rate by 
measuring the potential difference induced between a pair of electrodes by a flow of conducting liquid in the magnetic 
field. 3,4 This approach is now well developed and works reliably for common liquids like water, 5 but not so for liquid 
metals. Major problem in molten metals, especially at elevated temperatures, is the electrode corrosion and other 
interfacial effects, which can cause a spurious potential difference between the electrodes. 

The electrode problem is avoided by contactless eddy-current flowmeters, which determine the flow rate by sensing 
the flow-induced perturbation in an applied mag netic field. 6 - 7 The main problem with this type of flowmeters is the 
weak field perturbation which may be caused not only by the flow. We showed recently that the flow-induced phase 
shift of AC magnetic field is more reliable for flow rate measurements than the amplitude perturbation. 8 

Another contactless techniques for flow rate measurements in liquid metals is the so-called magnetic flywheel 
invented by Shercliff, 9 who prescribes a "plurality" of permanent magnets distributed equidistantly along the circum- 
ference of a disk, which is mounted on an axle and placed close to a tube carrying the liquid metal flow. The eddy 
currents induced by the flow across the magnetic field interact with the magnets by entraining them, which makes 
the disk rotate with a rate proportional to that of the flow. This type of flowmeter, described also in the textbook by 
Shercliff 2 and extensively used by Bucenieks, 10,11 was recently successfully reembodied under the name of the Lorentz 
force velocimetry (LFV). 12 

Recently, we suggested an alternative and much more compact design of such a flowmeter, which conversely to 
Shercliff 's flywheel uses just a single magnet mounted on the axle it can freely rotate around and magnetized perpen- 
dicularly to it. 14 We also introduced a basic mathematical model and presented first experimental implementation of 
this type of flowmeter. 15 When such a magnet is placed properly at a tube with the liquid metal flow, it starts to 
revolve similarly to Shercliff 's flywheel. But in contrast to the latter, which is driven by the electromagnetic force 
acting on separate magnets, the single magnet is set into rotation only by the torque. This driving torque is due 
to the eddy currents induced by the flow across the magnetic field. As the magnet starts to rotate, additional eddy 
currents are induced, which brake the rotation. An equilibrium rotation rate is attained when the braking torque 
balances the driving one, and this rate depends only on the flow velocity and the flowmeter arrangement, whereas 
it is independent of the electromagnetic torque itself. Thus, the equilibrium rotation rate is affected neither by the 
magnet strength nor by the electrical conductivity of the liquid metal provided that the friction on the magnet is 
negligible. This a major advantage of the single-magnet rotary flowmeter over the LFV approach, which relies on 
direct force measurements. 1 

In this paper, we present an extended theory of the single-magnet rotary flowmeter and compare it with exper- 
imental results. Two limiting cases of long and short magnets are analyzed using linear-dipole and single-dipolc 
approximations. We obtain simple analytic solutions for the force and torque on slowly moving and rotating mag- 
nets due to eddy currents in the layers of infinite horizontal extent and arbitrary depth. This allows us to find the 
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Figure 1: Schematic view of the single- magnet rotary flowmeter with the j/-axis directed out of the plane of figure. 

equilibrium rotation rate of the magnet at which the torques due to the translation and rotation balance each other. 
We also consider an active approach, where the force on the magnet is used to control its rotation rate so that the 
resulting force vanishes. This rotation rate, similarly to the equilibrium one, is proportional to the layer velocity and 
independent of its conductivity and the magnet strength. 

The torque on a magnetic dipole rotating about an axis normal to a thin sheet has been calculated by Smythe 
using an original receding image method. 16 Reitz uses this method to calculate the lift and drag forces on the coils 
of various geometries moving with constant velocity above a conducting thin plate. 17 The lift force, which at high 
speeds approaches the force between the coil and its image located directly below it, varies as the velocity squared in 
the low-speed limit considered in this paper. The drag force is found to vary inversely and directly with the velocity 
at low and high speeds, respectively. Palmer finds analytical expressions for the eddy current forces on a circular 
current loop moving with a constant velocity parallel to a thin conducting sheet. 18 The force on a rectangular coil 
moving above a conducting slab has been calculated numerically by Reitz and Davis using the Fourier transform 
method. 19 The same problem for the magnetic dipole of arbitrary orientation placed next to a thin slowly moving 
slab is addressed by Kirpo et al. 20 The force and torque on a transversely oriented dipole above a slowly moving plane 
layer of arbitrary thickness has been found analytically in the context of the LFV. 13 Fast computation of forces on 
moving magnets are of interested also for the eddy current force testing techniques. 21 

This paper is organized as follows. The following section presents two simple mathematical models of the single- 
magnet rotary flowmeter, which are used to calculate analytically the force and torque on the magnets moving and 
rotating slowly above a layer of infinite lateral extent. The limits of long and short magnets, which are approximated 
by linear and point dipoles, are considered in Sees. II B and II C, respectively. Section III presents the flowmeter 
implementation details and test results. The paper is concluded by a summary in Sec. IV. 



Consider a horizontally unbounded planar layer of electrical conductivity a occupying the lower half-space and 
moving as a solid body with a constant velocity Vq parallel to a permanent magnet placed at a distance h above 
its surface and rotating with a constant angular velocity f2 around an axis parallel to the surface and perpendicular 
to Vq. Velocities are assumed sufficiently low for the magnetic field of induced currents to be negligible compared to 
the field of the magnet. The origin of Cartesian coordinates is set at the center of the magnet with the x, y and z 
axis directed along Vq — ve x , ft — ile y , and downward normally to the surface, respectively, as shown in Fig. 1. In 
the following, two limiting cases will be considered in which the magnet will be assumed either much longer or much 
shorter than h. 



II. THEORY 



A. Formulation of problem 
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B. Linear-dipole model for a long magnet 

We start with a long cylinder magnetized perpendicularly to its axis about which it can freely rotate. 15 In this case, 
the magnetic field is approximated by that of a two-dimensional (linear) dipole with the vector potential A (r) = 
e y A (r) , which has only the y-component 

_ m ■ r _ hq m cos(-d - -dp) 

A \ r ) — 7; y~ ~ n ' y > 

Zir r z Zir r 

where /io = 4ir x 10~ 7 H/m is the vacuum permeability and m is the linear dipole moment, which is perpendicular 
to e y and directed at the angle do from the positive x axis; r is the radius vector from the magnet axis, and r = \r\ 
and $ are the cylindrical radius and the polar angle in the cylindrical coordinates around the y-axis. The magnetic 
field of linear dipole is 

B = V x A = -e y x VA (2) 

Eddy currents are induced by two effects: the translation of the layer and the temporal variation of the magnetic field 
due to its rotation. The latter vanishes in the co-rotating frame of reference, where the magnetic field is stationary, 
while the layer appears to move with the resulting velocity v = vq + v\, which contains also an apparent rotational 
motion of the layer v\ = fi x r opposite to that of the magnet. The density of eddy currents is given by Ohm's law 
for a moving medium 

3 = a(-V(p + v x B) (3) 

where ip is the electric potential. In this case, no electric potential is induced because the e.m.f., v x B, is both 
solenoidal and tangential to the surface. If the induced magnetic field is negligible as originally assumed, eddy 
currents can be represented as the superposition 

3 = Uo+.h)e y , 

where jo = —avo ■ VA and j\ = —ov\ ■ V^4 are the currents induced by the translation and rotation, respectively. 

It is important to note that in the approximation under consideration with a fixed magnetic field distribution, 
eddy currents are determined only by the instantaneous velocities. Besides that eddy currents are coplanar to the 
surface and, thus, mutually independent over the depth of the layer. Consequently, a layer of finite thickness may 
be represented as a semi-infinite one with zero velocity at z > hi, where hi is the distance of the lower boundary of 
finite-thickness layer from the magnet. This, in turn, is equivalent to the superposition of two semi-infinite layers with 
the second layer at z > hi moving oppositely to the first one at z > hi so that the resulting velocity vanishes at z > h 2 . 
In the following, this approach allows us to construct the solution for a finite-thickness layer by taking the difference 
of two half-space solutions, which is subsequently denoted by [X]^ , where X stands either for the force or the torque 
due to the eddy currents in a half-space. Moreover, by the same arguments, this approach can easily be extended to 
z-dependent velocity distributions v{z), for which general solution can be constructed as a superposition of solutions 
for thin layers moving with various velocities given by Jv(z)d z Xdz, where d z X = [X]% / [z]£ for h 2 -> h x = z and 
v = l. 

The linear force density experienced by an infinitely long magnet due to the layer translation, which according to the 
momentum conservation law is opposite to that acting upon the layer itself, can be written as — J s j x B ds = Foe Xl 
where the integral is taken over the xz-cross-section of the layer. The y-component of force is absent due to the 
y — > —y reflection symmetry. In the low-speed limit under consideration, when force varies linearly with the velocity 
there is also no z-component of force. This is the case because according to the linearity F — > —F when Vo — > —vo, 
while the latter transformation is equivalent to the rotation of the coordinate system by 180° around the z-axis, which 
leaves the z-components invariant. In polar coordinates with the surface defined by ri = h/sind, we obtain 

The linear torque density, which because of the aforementioned symmetries has only the y-component, can be found 
as — J s r x j'o x B ds = M e y , where 
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The linear force density due to the rotation, defined by — f s j\ x B ds = iqe x , is found as 

F 1= a f(e x -VA)( Vl -VA)ds= 



16nh 

The linear torque density due to the rotation, as that due to the translation above, has only the y-component 



Mi = a \ (r x VA)(vi ■ VA) ■ ds 



H 2 ) m 2 aVl 
4tt 2 



siiro; - ,) t) )d0 / y, 



(6) 



which is not defined for a semi-infinite layer because the last integral diverges. Nevertheless, expression (6) can still 
be evaluated for the layer of finite depth by substituting the infinite limit in the last integral by r-i = hi/ sin??, which 
results in 



8tt hi 1 ih2 



(7) 



The solution above becomes unbounded as hi — > oo and, thus, inapplicable to thick layers. This implies that for 
a half-space layer the induced magnetic field cannot be neglected however slow the rotation of magnet. The induced 
magnetic field is related to eddy currents by Ampere's law j = V x £?, which combined with expressions (2) and 
(3) leads to 

Mo<™i ■ VA + V 2 A = 0. 

As suggested by the external magnetic field (1), we search for the vector potential in the complex form 

A(r) = 3? lAiry^-^' . 

Then Eq. (8) for the complex amplitude A takes the form 

-ifiaaQA = r-^rA')' - r~ 2 A, 



(8) 



(9) 



(10) 



where the prime stands for the derivative with respect r. The general solution of Eq. (10) can be written as A(r) = 
CKi ((1 + i)r / d) , where C is an unknown constant to be determined by matching the induced and externally imposed 
magnetic fields; Ki(z) is the modified Bessel function of the second kind, 22 and d = y / 2/(/xoaO) is the skin depth due 
to the rotation. At distances r smaller than the skin depth (r/d <C 1), A& Cdj ((1 + i)r) is expected to approach the 
imposed field A = fi fh / (2irr) , which yields C = (1 + i)fiofh/(2ird) and 



A{r)= (l+^ Kl{{1 + i)r/d) 



(11) 



Substituting expressions (9) and (11) into integral (6), the torque on the magnet can be represented as Mi 



t£1 

2 JO 

after some algebra we obtain 



^ /; J(n(t?)) dtf, where I(r) = j r °° 3? 



A* A 



rdr and the asterisk denotes the complex conjugate. Using Eq. (10), 



d 2 

I( r ) = 
x 1 2 



'A* A'] = 5ft [K ((l + i)r/d)Ki((l - i)r/d){\ + i)r/d] . 

For the low rotation rates that satisfy r/S <C 1, we obtain 

!{r) « - O 2 5i [ln((l - i)r/d) + 7 ] = (£) ' (In(2d/r) - 7 ), 

where 7 = 0.577215 ... is Euler's constant. Finally, taking into account that ri($) = hi/ sin$ and J* ln(sin^) di? 
— 7rln2, we obtain 



Aigm'qfl j: 
8tt hi 



(12) 



where = de 7 is the effective skin depth, which owing to the low velocities under consideration is supposed to be 
large relative to hi. 
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Now we can use the results above to find the magnet rotation rate depending on the velocity of layer. First, 
if the magnet rotates steadily without a significant friction, the torques due to the translation and rotation are at 
equilibrium: [M + Mi]'^ 1 = 0, which yields 

°- 1 (£-£)"■£• < 13 > 

Note that this equilibrium rotation rate £1 depends neither on the magnet strength nor on the layer conductivity 
unless h 2 > d. In the latter case, the skin effect becomes important and, thus, /i 2 has to be substituted by d in the 
expression above. 

The velocity of layer can be determined also in another way by measuring the force on the magnet, as in the LFV, 
to control the magnet rotation rate so that the resulting force vanishes, [Fq + Fi] h * — 0. This results in 

n -3(sr + £)- ««) 

which again depends linearly on the layer velocity, but does not depend on its conductivity or the magnet strength. 



C. Single-dipole model for a short magnet 

In the other limiting case of a short magnet, when the distance to the surface is large or at least comparable to the 
size of magnet, the latter can be considered as a dipole with the scalar magnetic potential 

$(s;m) = -m-VG(x), (15) 

where m is the dipole moment and G(x) = (47r|a;|) _1 is the fundamental solution of Laplace's equation, which satisfies 
V 2 G = —S(x) with the Dirac delta function on the r.h.s. In the following, we use simplified notation e m • V = d m , 
where e m — m/m is the unit vector. Then the dipole magnetic field is given by B = — /io(V<I> — S(x)m), where the 
last term is added to ensure the solenoidality of B also at x = 0. The associated dipole current distribution is 

J = —V x B = V x 8m, (16) 
Mo 

which easily leads to the classical expressions for the force and torque used later on. In the following, we will be 
using also the spherical and cylindrical coordinates associated with the Cartesian ones in the usual way. As in the 
previous section, we change to the co-rotating frame of reference and consider the eddy currents due to translation 
and rotation separately. 



1. Translation 



For the translation with v n = ve x , the charge conservation V- j = applied to Eq. (3) results in Laplace's equation 
for <po 

VVo - 0. (17) 

At the surface z = h, the normal component of electric current vanishes: 

d z f = vB y . (18) 

In order to find the induced electric potential, firstly, it is important to notice that B being a free-space magnetic 
field satisfies the Laplace equation itself. Consequently, the Cartesian components of B satisfy this equation, too. 
Secondly, if B y in BC (18) satisfies Laplace's equation, then 



v J B y dz = Li a mv J dy m Gdz. (19) 



satisfies not only BC (18) but also Eq. (17) because the integration along a straight line, similarly to the differentiation, 
are interchangeable with the Laplacian. By the same argument, we can interchange integration and differentiation in 
expression (19), which yields 

(po = fi mvdy m H, (20) 
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where H = jGdz = j^(Q (cos9) + q(r)) and Qq(z) = lnyi±| i s the zeroth degree associated Legendre function 

of the second kind; 22 8 is the spherical polar angle from the positive £-axis (see Fig. 1) and r is the corresponding 
cylindrical radius. The "constant" of integration q(r), which similarly to the first term Qq(cos9) is supposed to be 
axisymmetric and also to satisfy the Laplace equation, is chosen to regularize H at r — > by removing the logarithmic 
singularity Q a = In cot | <~ — lnr for z > 0. This results in q(r) = lnr and 

H=±-(Q (z/R) + lnr), (21) 

where R = \x\ is the spherical radius. For a transversal dipole (e m — e z ), considered also by Thess at a/. 13 , solution 
(20) reduces to 

fa = HovmdyG = —iiov<&{x]me y ), (22) 

where the last term represents the magnetic potential of a dipole aligned with the y-axis. For a general dipole 
orientation, expression (21) substituted into solution (20) after some algebra yields 



^ = 17 



(R — z z \ 2m z rsva.d — m v z 
VR^ ~2R3j( mx Sm ~ my C ° S ' + _ 



2i? 3 



(23) 



where $ is the azimuthal angle from the positive a;-axis in the xy-plane. 

The electric potential distribution (23) allows us to calculate the force acting upon the magnet, which, as noted 
above, is opposite to that acting upon the layer, i.e., F = — J v f dV, where the integral of the Lorentz force density 
/ = j x B is taken over the layer volume V. For the longitudinal force component, we obtain 

F , x - ^ m 2 avJ s [(d 2 ym H) 2 - (d m G) 2 ] + J\d 2 xm G) 2 dzj ds 

^^ml + ml + Aml), (24) 



512tt/i 3 

where the first integral is taken over the surface S at z = h and can be swapped with the second one over the layer 
depth. For a transversal dipole (m x = m y = 0), the expression above coincides with that of Thess et at 13, as well as 
with the result of Reitz 17 for a thin sheet in the limit of a slowly moving dipole. As seen, the force is the strongest 
on a transversal dipole and reduces on longitudinal and spanwise dipoles by factors of | and \, respectively. Thus, 
the force (24), in contrast to force (4) for a long magnet, varies with the dipole orientation in the xz-plane. On a 
horizontally inclined dipole (m x m y ^ 0), there is also a spanwise force component 

F , y = ^m 2 ^|J(a 2 m if)(a 2 y G) + ^(aLG)(a 2 m G)dz|d s 



\i\ov 



256^™- 

But there is no vertical force whatever the dipole orientation: 



(25) 



F^ z = n 2 av ^ [ (d 2 m G)(d 2 zm G) dsdz = 0. (26) 
Jh Js 

This is consistent with the result of Reitz 17 stating that in the low-speed limit the lift force is proportional to the 
velocity squared, while the approximation under consideration takes into account only the part of force proportional 
to the velocity. 

Alternatively, the force and torque acting on dipole can be found using the associated current distribution (16) and 
the induced magnetic field b, which lead straightforwardly to the classical expressions 23 

F= [jxbd 3 x = (m-V)b, (27) 
Jv 

M = / x x J x bd 3 x = m x b, (28) 
Jv 

where the integrals are taken over the space V above the layer. This is a bit longer but algebraically more straight- 
forward approach, which will be pursued in the following. 
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In order to find the induced magnetic field b, it is important to notice that solution (19) satisfies condition (18) not 
only at z = h but at any z. Thus, the ^-component of the current is absent not only at the surface but throughout 
the whole layer. Then V • j = and j z = imply 

j = V x xjje z , (29) 

where ip is the electric stream function, whose isolines coincide with the eddy current lines. Substituting this expression 
into Eq. (3) and taking the z-component of the curl of the resulting equation, we obtain 

(V 2 -d 2 z )^ = -cre z -Vx(vxB). (30) 

Note that in contrast to the electric potential in Eq. (17), no boundary conditions are required for ip at the surface. 
This is because, firstly, Eq. (30) contains no derivatives in z and, secondly, the absence of the z-component of current 
at the surface is explicitly ensured by expression (29). 

For translational motion with v n = ve Xl the r.h.s of Eq. (30) takes the form avd x B z = —[i <Tvd xz $ 7 which by the 
same arguments as above satisfies Laplace's equation and equals to —d z ipo with 

V>o = -n mavdl m H, (31) 

which satisfies the same equation and, thus, represents the solution of Eq. (30). For a transversal dipole (e TO = e z ), 
the general solution (31) simplifies to 

?Ao = —HomavdxG = iiQ(jv<f>{x; me x ), (32) 

where the last term represents the magnetic potential of a longitudinal dipole. Thus, in this case, the eddy current 
lines coincide with the isolines of electric potential (23) rotated by 90° about the z-axis. For a longitudinal dipole 
(e TO = e x ), the comparison of solutions (31) and (20) shows that the eddy current lines and the electric potential 
isolines are swapped with the corresponding distributions induced by a spanwise dipole (e m = e y ) and rotated by 90° 
about the z-axis. For an arbitrarily oriented dipole, expression (21) substituted into the general solution (31) after 
some algebra yields 



R-z z \ . 2m z r cos - m x z 
~lh^~ W?) (m x cos 2^ + 7% sin 2tf) + ^ 



(33) 



Figure 2 shows the isolines of the electric potential (23) and those of tpo, which represent the eddy current lines, 
in the xy-plane [z = const) for three basic dipole orientations along the x-, y- and z-axis. Note that the patterns 
are self-similar in the xy-plane with the characteristic length scaling directly with z. Namely, these distributions are 
functions of spatial angles only, while according to (23) and (33) the magnitude of both ip and f/'o falls off as ~ z~ 2 . 

In order to satisfy the solenoidality condition V • b = 0, the induced magnetic field is sought as b — V x a, where 
a is the vector potential, which is supposed to satisfy the Coulomb gauge V • a = 0. Then Ampere's law leads to 
V 2 a = —Hoj, which, in turn, results in a(x) = J v j(x')G(x — x') d 3 x'. Substituting expression (29) into the last 
integral, after some algebra we obtain a = [i V x x e z, where 



[ il>(x')G(x-x')d 3 x' 
Jv 



is the same as used by Thess et al. 13 The previous expression implies that a, similarly to its source j, has no 
z-component. Then the induced magnetic field can be written as 

b = /i V x V x X e z = Vo(e z ip - V</>), 

where = — d z \ represents the scalar magnetic potential, which completely defines b outside the layer, where 
%jj = 0. Further, using the tensor notation with the Einstein summation convention, expressions (27) and (28) can be 
represented in terms of \ as 

F t = fJ, mjX;ij3, (34) 
Mi = ^ijk l m j x ;k3 , (35) 

where Cijk is the anti-symmetric tensor and the subscript after the semicolon denotes the differential with respect 
to the corresponding coordinate with indices 1, 2 and 3 standing for the x,y and z, directions, respectively. Taking 




Figure 2: 
represent 
z-axis (c) 



Isolines of the electric potential ipo (dashed) and of the stream function tpo (solid) in the xy-plane, where the latter 
the eddy current lines, induced by translation in the magnetic field of the dipole aligned with the x- (a), y- (b), and 



into account the symmetry of G(x — x') with respect to the interchange of observation and integration points, the 
derivatives of x at the dipole location x = in expressions (34) and (35) are found as 



X-ij — 
X'.ijk — 



\ ^(x)G.. ij (x)d 3 x, 
Jv 

ip(x)G. ijk (x)d 3 x, 



(36) 
(37) 
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where 



■ij{x) = ■ 



G-ij k {x) = - 



4tt|x| 5 

47r|a;| 7 



Integrals (36) and (37) can be evaluated analytically using, for example, the computer algebra system Mathematica, 24 
which also allows us to carry out all other analytical transformations. In such a way, we firstly verify that Eq. (34) 
with expression (31) substituted into integral (37) indeed reproduces previous results (24), (25) and (26). Secondly, 
Eq. (35) with the same expression for tp substituted in integral (36) results in 

2 

M ° = 128ttI 2 (^ mxm v ex + ( m * + m l) e v - m y m z e z ) , (38) 



which for a transversal dipole again coincides with the results of Thess et al. 
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2. Rotation 



For a solid-body rotation with the velocity V\ = — f2 x x, which appears in the co-rotating frame of reference, the 
charge conservation V • j = applied to Eq. (3) results in 



VVi = -2QB y , (39) 
while the vanishing of the normal current component at the surface z = h requires 

d z (f! = -ClzB y . (40) 
In order to solve this problem, firstly, it is important to notice that since J B y dz satisfies the Laplace equation, 

if i = flz J By dz = ^mVLzdy m H (41) 

is a particular solution to Eq. (39). Comparing expressions (41) and (23) shows that (pi = (£lz/v)ipo, where ifo is 
given by solution (20). Searching for the solution as f\ = <fi + ifi, reduces Eq. (39) to the Laplace equation for <fi, 
while BC (40) takes the form 

d z fi = -fiQmfldy m H. 

By the usual arguments, we obtain 

<fi = -fiomfldymHi, (42) 

where Hi = J H dz = (RQi(cos9) + z\nr) = zH — R 2 G and Qi(z) = zQo(z) — 1 is the first order associated 
Legendre function of the second kind. 22 Using the expressions above, after some transformations we obtain 



f i = —^o^ ( e x ■ m x VH + ymd m G) 



47T 



m z (R? — z 3 ) sin i? (m x sin2i9 — m y cos2i9)r 2 — m y (R 2 + z 2 ) 
Wr + 2^3 



(43) 



which is the electric potential in the co-rotating frame of reference. 

As before, condition (40) is satisfied not only at z — h but any z and, thus, the z-component of the current vanishes 
throughout the layer. Consequently, the electric current can again be expressed as (29) with ip\, for which Eq. (30) 
takes the form 

(V 2 -d^^aniBx + eyrxVBz). (44) 
Since the r.h.s. of Eq. (44) satisfies the Laplace equation, after some algebra we obtain 



-an (^jj B x dz 2 + e y -rxV JjB z dz 2 -2d x JJJ B z dz 3 ^j 



^ tinuO f m z (R — z) cos — m x r^, 
/i aQe y -mxV5 =^^(— ^ — )■ (45) 
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Figure 3: Instantaneous isolines of the electric potential y>i (dashed) and of the stream function ipi (solid), which represent 
the eddy current lines, induced by a solid-body rotation in the magnetic field of the dipole aligned with the x- (a), y- (b), and 
2-axis (c). 



The eddy current lines in the xy-p\ane, which are the isolines of ipi, along with the isolines of the electric potential 
(43) are shown in Fig. 3. Again, the patterns are self-similar and scale directly with the distance z from the dipole, 
while the magnitude of both and ipi falls off as <~ z^ 1 . In contrast to the translation considered above, the eddy 
currents induced by rotation are time-dependent and vary periodically as the dipole orientation changes from the 
x- to the z-axis, which are shown in Figs. 3(a) and (c), respectively. Only the electric potential, but no current, 
is induced by a solid-body rotation when the dipole is aligned with the axis of rotation (y). In his case, when the 
magnetic field is symmetric about the axis of rotation, the induced e.m.f. is irrotational and, thus, compensated by 
the electric potential gradient, which equals to 

—vi x B = -Qpe x V x (Ae) = -SlV(pA), (46) 
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Figure 4: Instantaneous rotation rate (51) versus time (solid) and orientation (dashed) for a negligible inertia of the magnet. 



where p = \ fl x x\/il is the cylindrical radius from the axis of rotation, e = f2xa;/|nxx|is the azimuthal unity 
vector, and A is the azimuthal component of the vector potential, which can be used to describe a general axially 
symmetric poloidal magnetic field. For the dipole field, we have 

A(x) = p [ J(x')G(x - as') dV = m x VG(x) = A(x)e, 
Jv 

where J is the associated dipole current (16) and A(x) = —\m x x\G(x)/\x\ 2 . Substituting this into expression (46) 
and equating it to the gradient of the electric potential ipi, we obtain 

^m y n(R 2 - y 2 ) 

^-~ Q P A ~ > 

which is the free-space electric potential in the co- rotating frame of reference shown in Fig. 3(b) and coinciding with 
(43) when m x = m z = 0. This potential vanishes in the laboratory frame of reference, where the magnetic field is 
invariant with respect to the rotation around the symmetry axis. 

Further, using expression (45) for ipi in integrals (36) and (37), which can be evaluated together with Eqs. (34) 
and (35) in the same way as for the translation in the previous section, we obtain 

Fi = (~( m x + ™%)e x ~ m x m y e y + m x m z e z ), (47) 

M ^ = (2m x m y e x - (2m 2 x + m 2 z )e y + m y m z e z ). (48) 

647m 

The first point to note is that both the force and torque vanishes when dipole is aligned with the axis of rotation 
(m x = m z = 0), which, as discussed above, is due to the absence of eddy currents in this case. Secondly, the 
longitudinal (x) component of force (47), in contrast to that due to the translation given by expression (24), is 
independent of the dipole orientation in the xz-plane and depends only on the magnitude of the dipole moment in 
this plane: m 2 = m 2 +ml. Thirdly, there is also a non-zero transversal (z) force, which in contrast to the longitudinal 
one is purely oscillatory and varies periodically with the dipole orientation in the xz-plane as m x m z — ^rh 2 sin2-# , 
where $o is the poloidal angle of the dipole orientation in the xz-plane from the positive x-axis, which is the same 
as that for the linear dipole in expression (1) (see Fig. 1). Moreover, when dipole is inclined to the axis of rotation 
(m x m y ^ 0), also a spanwise (y) force appears, which similarly to the transversal one alternates periodically around 
zero as m x m y — mm y cos z9o with the dipole rotation. As seen from expression (48), inclination also gives rise to 
alternating torque components around both the x- and z-axis. The y-component of torque, which is negative and, 
thus, opposing the rotation, has not only a constant but also an alternating part, which varies with the double 
frequency of rotation as —(2ml + m l) = — f 'ti 2 (l — 3 cos2$ ) • 

Thus, in contrast to the long magnet in §11 B, there is no steady rotation rate, which could balance the y-component 
of the constant torque (38) due to the translation. The variation of the magnet rotation rate in response to oscillatory 
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torque is constrained by the inertia of magnet. For the oscillatory part of the rotation rate 17 to be negligible compared 
to the mean one £7, we need 

~ - — ~ ™ < 1 (49) 

where / is the moment of inertia of magnet. If this condition is satisfied, the oscillatory component may be neglected 
in the balance of torques (38) and (48) around the axis of rotation, which is the y-axis. For a layer of finite thickness, 
this results in 

= U±- + f). (50) 



3 V hi h 



The expression above differs only by a factor of ^ from result (14) for a long magnet in the case of a vanishing 
longitudinal force. In the opposite limit of a negligible inertia, the torque balance yields 

n = n /(i - |cos2?? ). (51) 

This instantaneous rotation rate £1 — dt'&o, which is shown in Fig. 4 versus both the time t and orientation i9 , 
alternates between |f2 and |Qo of the mean value (50) Cl = i?oAlt_>.oo • The last relation follows from the integration 
of Eq. (51) as i?o — ^ sin2#o — which also defines SI parametrically versus t in Fig. 4. Note that the oscillations 
of the rotation rate do not affect its mean value il , which is defined by (50) independently of condition (49). 

Alternatively, when the magnet rotation rate is actively controlled so that to balance the longitudinal force (24) 
with the ^-component of force (47), the instantaneous rotation rate is 

fl = d t $ Q = — U7- cos2tf ), (52) 



where the temporal mean value 



2 \hi h 2 hi + h 2 
follows from the solution of differential equation (52), which can be written as 

tan(tf (*)) = ^j-tan(fiii). 
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III. IMPLEMENTATION AND TESTING 



The laboratory model of the single-magnet flowmeter shown in Fig. 5 consists of a cylindrical SmCo-type permanent 
magnet with diameter 2R = 24 mm and height L = 35 mm, magnetized perpendicularly to its axis with 0.5 T surface 
induction, which holds up to 230°C. The magnet is mounted on a stainless steel axle held by a housing made of 
"MACOR" ceramic, which ensures a low mechanical friction and can withstand up to 800°C. 

The flowmeter was tested on a sodium loop using about 90 liters of molten Na at 220° C with the electrical 
conductivity of 9 x 10 6 Sm _1 . The flow was driven by a linear electromagnetic pump, which provided a maximum 
velocity of 1.5 ms' 1 at the flow rate of 31s -1 in a stainless steel duct of 45 x 45 mm 2 cross-section. 

The frequency of magnet rotation was determined with an inductive magnetic proximity sensor (SICK MM12- 
60APS-ZU0), which was fixed at the distance of 50 mm perpendicularly to the side surface of magnet at its midheight 
as shown in Fig. 5. The sensor briefly switched off and then on again as the magnetic field along the axis of the 
sensor changed its direction. The frequency of sensor transition from the off to the on state, which happens twice 
per revolution, was measured with Keithley 2000 Multimeter using the reciprocal frequency counting techniques with 
the measurement gate time set to 1 s. The frequency was also monitored with Tektronix TDS 210 oscilloscope. After 
changing the pump current, the flow was allowed to develop for 2 min, and then six measurements were taken with 
the intervals of approximately 30 s. The standard deviation in the measured rotation frequencies was typically a few 
tenths of Hz. 

The rotation rates measured at two gap widths between the magnet and the duct are shown in Fig. 6(a) depending 
on the average sodium velocity, which was determined using Ultrasonic Doppler Velocimetry 25 with the resolution of 
1 cms" 1 . 26 Although the rotation rate is seen to increase as the magnet is approached to the duct and to vary nearly 
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Figure 5: Single-magnet rotary flowmeter with a SmCo-type magnet (1) of 24 mm and 35 mm in diameter and height, respec- 
tively, installed at the side of a rectangular stainless steel duct (2) of 45 x 45 mm 2 in cross-section carrying the liquid sodium 
flow. The rotation rate is measured with an inductive magnetic proximity sensor (3) and also monitored with a photoelectric 
barrier sensor (4). 

linearly with the flow velocity in agreement with the theory, the linear fit shows a certain zero offset. Obviously, there 
is some opposing force, which has to be overcome for the rotation to start. This requires the flow velocity of at least 
0.2 ms -1 . First, such an opposing force is caused by the dry friction in bearings. Second, the magnet is oriented by 
the Earth's magnetic field, which has a stronger effect than the friction and also has to be overcome for the rotation to 
start. Third, even stronger effect is caused by the stainless steel walls of the duct, which are weakly ferromagnetic. 27 
The magnet, when approached to the duct, was observed to turn with the dipole moment perpendicular to the wall. 
It is also important to note that the magnet, when turned by hand away from the equilibrium orientation, returned 
to it monotonically without oscillations. This implies the inertia of magnet to be small relative to the electromagnetic 
drag force. Thus, the inertia cannot significantly contribute to the overcoming of the orienting torque due to the 
wall magnetization. However, when magnet rotates slowly, the orienting torque causes the rotation rate to fluctuate, 
which shows up as an increased scatter in the measured values seen in Fig. 6 at low flow velocities. 

When the zero offset is removed and the remaining rotation rate is multiplied by the distance h = R + d w + d 
between the magnet axis and the liquid metal, where R = 12 mm is the magnet radius and d w = 3 mm is the thickness 
of stainless steel wall, we obtain the relative rotation rate k = Qh/v, which represents a rescaled slope coefficient for 
Fig. 6(a) and is plotted in Fig. 6(b) versus the velocity. The short-magnet solution (50) in two limiting cases of a 
semi- infinite (/12 = 00) and a thin [ho, = hi) layer yields k = g and k = |, respectively, which are smaller than the 
measured values k sa 1 seen in Fig. 6(b). On the other hand, the long-magnet solution (13), which appears more 
adequate for the experimental setup, yields k < i for hi < h^. Such quantitative differences between the experiment 
and theory are not surprising given the simplifications underlying the latter. First, theoretical model does not take 
into account the finite width of the duct. But this alone could hardly explain the high rotation rate of the magnet 
at which its field travels faster than the layer. This apparently being the case in the experiment implies the presence 
of significant velocity gradients, which are also ignored by the theory, but could be taken into account as outlined 
in Sec. II B. Note that strong velocity gradients can be caused by the magnet itself, which due to its large size and 
strength may act as a magnetic obstacle partially blocking the flow and so increasing its local velocity. 

IV. SUMMARY AND CONCLUSIONS 

We have presented a theory of single-magnet rotary flowmeter for two limiting cases of long and short magnets, 
which were modeled as linear and single dipoles, respectively. Simple analytical solutions were obtained for the force 
and torque on slowly translating and rotating magnets due to eddy currents in layers of arbitrary depth and infinite 
horizontal extent. The velocity was assumed to be constant and the motion so slow that the induced magnetic field 
could be neglected. The latter assumption was not applicable to the long magnet rotating above a conducting half- 
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Figure 6: Measured (a) and rescaled (b) magnet rotation rates along with the best linear fits versus the average velocity for 
d — 8 mm and 11mm gaps between the magnet and the duct with the distance between the magnet axis and liquid metal 
h — R + d w + d, where R = 12 mm is the magnet radius and d w = 3 mm is the thickness of stainless steel wall. 



space. In this case, to obtain a finite braking torque, the skin effect due to the induced magnetic field had to be taken 
into account however slow the rotation. For a single dipole of arbitrary orientation, compact analytical solutions were 
obtained in terms of both the electric potential and stream function induced by the layer translation and rotation 
in the co-rotating frame of reference. The electric stream function was used further to find the scalar potential of 
induced magnetic field at the dipole location, which resulted in simple expressions for the force and torque on the 
dipole. Eventually, we found the equilibrium rotation rate at which the driving torque due to the layer translation is 
balanced by the braking torque due to the magnet rotation. An alternative approach was also considered, where the 
force on the magnet could be used to control its rotation rate so that the resulting force vanishes. In either case, the 
resulting rotation rate is directly proportional to the layer velocity and inversely proportional to the distance between 
the magnet and the liquid metal. These results were found in a qualitative agreement with the measurements on 
the liquid sodium flow. A more accurate quantitative agreement with experiment is limited due to the substantial 
approximations underlying the theoretical model, which neglects the finite lateral extension of the layer as well as the 
spatial and temporal variations of the velocity distribution. 

In conclusion, note that the resulting rotation rate is independent of the magnet strength and the electrical con- 
ductivity of the liquid metal provided that the mechanical friction or other external effects are negligible compared 
to the driving torque. This is the main advantage of rotary flowmeter over the LFV. 12 
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